The Low-Lying Dirac Spectrum of Staggered Quarks 
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We investigate and clarify the role of topology and the issues surrounding the e-regime for stag- 
gered quarks. We study unimproved and improved staggered quark Dirac operators on quenched 
lattice QCD gluon backgrounds generated using a Symanzik-improved gluon action. For the im- 
proved Dirac operators we find a clear separation of the spectrum into would-be zero modes and 
others. The number of would-be zero modes depends on the topological charge as predicted by 
the continuum Index Theorem, and the expectation values of their chirality are large for the most 
improved actions (~ 0.7). The remaining modes have low chirality and show clear signs of cluster- 
ing into quartets that become degenerate in the continuum limit. We demonstrate that the lattice 
spacing and volume dependence of the eigenvalues follow expectations. Furthermore, the non-zero 
modes follow the random matrix theory predictions for all topological charge sectors. The values 
of the chiral condensate extracted from fits to the theoretical distributions are consistent with each 
other, and with the results obtained from the total density of eigenvalues using the Banks-Casher 
relation. We conclude that staggered quarks respond correctly to QCD topology when both fermion 
and gauge actions are improved. 

PACS numbers: ll.15.Ha, 12.38.Gc 



I. INTRODUCTION 

The low energy regime of Quantum Chromodynam- 
ics (QCD) exhibits a rich and interesting phenomenol- 
ogy, including the 1/4(1) axial anomaly, chiral symmetry 
breaking and the topological properties of the theory. A 
crucial step in elucidating these effects is understanding 
the low-lying eigenvalue spectrum of the Dirac opera- 
tor. There are a number of detailed predictions of the 
properties of these low-lying modes, such as the exis- 
tence of an Index Theorem, the Banks-Casher relation, 
and the distribution of the first few eigenvalues in fixed 
topological charge sectors. Such effects are, however, in- 
herently non-perturbative and can only be studied fully 
using techniques like lattice Monte Carlo simulation. 

Universality decrees that any correct discretization of 
QCD must be correct close enough to the continuum 
limit. If we want to make phenomenological predictions, 
however, it is crucial that the QCD lattice Dirac opera- 
tor also has the correct low-lying spectrum at the lattice 
spacings typically simulated. By this we mean that the 
deformations to the eigenvalue spectrum brought about 
by the discretization must be small. In this case the rel- 
evant scale for comparison is the light sea quark mass, 
TO„,d = 0(5 MeV). 

In this paper we show that these deformations are 
small enough for one particular class of lattice fermions 
that are widely used in large scale lattice simulations: 
improved staggered quarks. Working at similar lattice 
spacings in the quenched theory (i.e. pure Yang-Mills 
gluodynamics) we show that the low-lying QCD Dirac 
eigenvalue spectrum very closely reproduces the contin- 
uum features. In particular, we show that improved stag- 
gered fermions do respond correctly to the gluonic topo- 



logical charge, and discuss why some confusion on this 
issue exists in the literature. 

To reach these conclusions we test a number of pre- 
dictions, quantitative and qualitative, for the low-lying 
modes: the Atiyah-Singer Index Theorem, that relates 
the zero modes to the topological charge of the gauge 
field; the Banks-Casher relation, linking the chiral con- 
densate to the density of eigenvalues at the origin; the 
universality of the low-lying eigenvalue spectrum in the 
£-regime, and its connection to Random Matrix Theory. 
In addition, we change both the lattice spacing and vol- 
ume and show that the variation of the spectrum matches 
theoretical expectations. 

This paper expands and refines the work in P, ■ We 
include much more extensive data, both on larger vol- 
umes and finer lattices. Related work has been presented 
in 

The structure of the paper is as follows. In Sec. ITTlwe 
describe our expectations for the Dirac eigenvalue spec- 
trum, both in the continuum and with lattice staggered 
fermions. We also review the existing literature. We 
describe our methodology in Sec. IIIII and discuss the 
features of the spectrum in Sec. II VI In Sec. we com- 
pare our results for the low-lying eigenvalues with the 
universal predictions from random matrix theory. We 
summarise and conclude our study in Sec. IVII 



II. THEORETICAL BACKGROUND 

In this section we review the theoretical expectations 
for the spectrum of the Dirac operator in the continuum 
and lattice staggered cases. 
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A. Continuum QCD 

A single, continuum fcrmion species is described by the 
massless, gauge covariant Dirac operator, which is anti- 
Hermitian with a purely imaginary eigenvalue spectrum: 



Pfs = tXsfs , A, e M . (1) 
Choosing orthonormalised eigenvectors, satisfying 

/] /* = Sst , (2) 
the chirality of an eigenvector is defined by 

Xs = fhbfs ■ (3) 
The Dirac operator also anticommutes with 75, 

{^,75} = 0, (4) 

implying that the spectrum is symmetric about zero: 

sp(^) = {±i\s, A, e M} . (5) 

If As ^ 0, then 75/s is also an eigenvector with eigenvalue 
—iXs, and zero chirality, Xs = 0. The zero modes, As = 0, 
can be chosen with definite chirality, Xs = if- 

1. Topological charge and the Index Theorem 

Smooth continuum gauge field configurations can be 
classified by their winding number, which is an integer 
given by the expression: 



Q = [ d'^x e^^ar Tr Ff_,i,{x)F„r{x) 

AZtt^ J 



(6) 



If we now consider the Dirac operator in the given gauge 
background, and we denote by n± the number of zero 
modes with a given chirality, their difference is related 
to the to;^logical charge via the Atiyah-Singer Index 
Theorem 



Q = m Tr 



75 



P{A) + TO 



= ^Xs=n+- n- 



(7) 



where to is the quark mass (we keep the sign conventions 
ofRef. il,]). 



2. The chiral condensate and the Banks-Casher relation 

The bare chiral condensate in a given gauge field back- 
ground, A, is given by the trace over the quark propaga- 
tor: 



1 



(V'V'M) = -T7 Tr 



1 



V p{A) + TO 



(8) 



In the limit of to ^ it acts as an order parameter for 
the spontaneous breaking of chiral symmetry. Expanding 
in eigenmodes of If) and using the ±iAs symmetry. 



V 



(9) 



Averaging over a full ensemble of gauge fields, the sum 
over s can be replaced by an integral over eigenvalues 
with measure given by the spectral density: 



P(A)^:^^^(A-As 



(10) 



The integral has a smooth limit as to —s- 0: 

S = - l~4)i)(m = 0)) = limm [ dX ^f^'^'' , 



^p(O) , 



(11) 



which is the Banks-Casher relation. We expect the chi- 
ral condensate E to be non-zero when the chiral sym- 
metry is spontaneously broken. This requires that the 
large volume limit is taken before to ^ to ensure that 
the discrete sum over eigenmodes can be replaced by the 
smooth integral. 



3. The e-regime and spectral universality 

The e-regime of QCD occurs when the theory is reg- 
ularised in a finite volume such that the typical system 
size L satisfies 



e-regime: 



(Aqcd) ^ < L < {M^y 



(12) 



In such systems, the volume is large enough that we may 
derive an effective chiral field theory for pions by integrat- 
ing out all other hadronic states. The typical Compton 
wavelength of the particles so removed is (Aqcd)^^- At 
the same time, the system is very small compared to the 
Compton wavelength of the pions, so we can further inte- 
grate out all their kinetic modes to leave us with a theory 
describing just static pions. 

It can then be rigorously shown that the remaining 
finite volume partition function exactly corresponds to 
that of a random matrix theory (RMT) with the same 
chiral symmetries as QCD 0, This universality 

0, allows us to predict the distributions of low- lying 
non-zero eigenvalues of the Dirac operator using RMT 
[Til ,i12 ] . The relevant universality class is determined by 
the chiral symmetries of the finite volume partition func- 
tion, and hence of QCD; for fermions in the fundamental 
representation of the gauge group with Nc > 3 the appro- 
priate class is the chiral Unitary Ensemble (chUE) with 
Dyson index /3 = 2. 

Rather than the absolute eigenvalue spacings, RMT 
predicts only the relative spacings of the eigenvalues in 
the form of microscopic (or unfolded) spectral densities 
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Ps{z). These detail the universal distribution of one or 
more eigenvalues. There are separate predictions for each 
sector of topological charge. The absolute spectral den- 
sities include a QCD-specific scale factor based on the 
chiral condensate: 



(13) 



Comparison of the Dirac spectrum with RMT therefore 
provides a method of determining the chiral symmetries 
of QCD, of probing the role of topology, and of measuring 
the chiral condensate. 

For the combined spectral density of all eigenvalues in 
a given topological charge sector RMT predicts 



J.+iiz)J.-i{z)] , (14) 



where Ju{z) are Bessel functions and = \Q\. As these 
and all other predictions depend only on \ Q\, we improve 
the statistical accuracy of the comparison with the Dirac 
spectrum by combining the ±|(5| sub-ensembles. 

There are also separate predictions for the fc*^ eigen- 
value [Hill, e.g. k = l: 

§e"^V4 , ^ 0) , 

PiW- <( §/2(^)e-^'/\ (^^ = 1), 

|(/|(z)-/3(z)/i(z))e-^/\ (z.^2). 



Expressions for higher k are given in Ref . ^3 • 



(15) 



B. Lattice staggered fermions 

The massless, gauge-invariant, One-link staggered 
Dirac operator on a d = 4 dimensional Euclidean lattice 
with spacing a is 



1 



2auQ 



r]^{x) [U^{x)Sx+fL,v - H.c] (16) 



(17) 



with uq an optional tadpole-improvement factor given, in 
our case, by the fourth root of the mean plaquette. Ip is 
antihermitian and obeys a remnant of the continuum 75 
anticommutation relation: 

{p, e} = , with e{x) = (-1)^^=1 . (18) 

As in the continuum case, its spectrum is therefore purely 
imaginary, with eigenvalues occurring in complex conju- 
gate pairs: 



sp(^) = {±a„A, e 



(19) 



In fact, if fs is an eigenvector with eigenvalue iXg, then 
e fs is also an eigenvector with eigenvalue —iXg- 

The corresponding action describes Nt = 2'^^^ = 4 
"tastes" of fermions which interact via unphysical "taste 



breaking" interactions that vanish in the continuum limit 
as a^. In the limit a there is an SU{Nt) SU{Nt) 
chiral symmetry, and the spectrum is therefore an exact 
A^4-fold copy of that described in Sec. Ill Al 

At finite lattice spacing the chiral symmetry group is 
reduced to L/(l) C/(l) and we do not expect to see this 
picture, for example there will not be an exact Index 
Theorem anymore. The relevant operator for the Index 
Theorem must be a taste-singlet, 75^, to couple to the 
vacuum correctly [T5|. It is a point split operator in all 
d-dimensions, which can be made gauge invariant by the 
insertion of the appropriate U fields. It does not anticom- 
mute with the staggered Dirac operator, {_p,75*} ^ 0. 
This also implies that the chirality of the eigenmodes 
will not be exactly ±1 or 0. 

Exact zero eigenvalues only occur accidentally (and are 
completely excluded at non-zero fermion mass), and all 
the eigenmodes will in principle contribute in Eq. Q |l6j| . 
Due to the simple properties 



0, 



(20) 
(21) 



each pair of eigenvectors {fs, e fs} has the same chirality 

Xs- 

(efsiiri^fs) - {f.\e^ir^\fs) = (fsbrifs) ■ m 

All of the preceding considerations depend only on sim- 
ple symmetry properties, and apply equally well to the 
improved staggered operators we have used in our work. 

Early studies of the One-link staggered operator and 
unimproved, Wilson g aug e action for coarse, thermalized 
lattices were done in 0, 0, 0, 0, |23| . There was no 
clear separation in the low-lying modes, either in eigen- 
value or chirality. The low-lying One-link eigenvalues 
also did not follow the universal predictions for contin- 
uum QCD. In particular they showed no dependence on 
the topological charge sector '2l'. "2?. '23'. "2112511 . 

In this study, however, we are interested in whether im- 
provement of the gauge action and Dirac operator allow 
the continuum-like spectrum to be seen. 

Certainly, some continuum features emerge using very 
smooth gauge backgrounds (obtained by either repeated 
cooling or discretising classical instanton solutions) 0, 
IT^ . By contrast, here we are interested in thermalized, 
unsmoothcd configurations at lattice spacings (and gauge 
actions) typical of lattice simulations, i.e. a « 0.1 fm. We 
conjecture that these failures of continuum predictions 
follow from the use of coarse lattices and unimproved ac- 
tions, resulting in large taste changing interactions that 
ruin the continuum chiral symmetries. We shall show 
that the improved staggered operators introduced be- 
low allow all the continuum features of the spectrum to 
emerge. 

That improvement allows this to happen at such lattice 
spacings is not entirely unexpected. QCD predicts that 
the topological susceptibility is suppressed as we reduce 
the sea quark mass. This is seen in simulations using 
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improved gauge and staggered fermion actions 

With One-link staggered sea quarks no such 
variation is seen 1^ • Other evidence for 

the beneficial effect of improved staggered operators was 
seen in the Schwinger model .3^ 3It> i^j i as well 
as in QCD [IIIIMEI- 



1. Improved staggered fermions 

The continuum chiral symmetry is broken by the (un- 
physical) taste-changing interactions. At leading order 
these involve exchange of an ultraviolet gluon of momen- 
tum q n/a. Such interactions are perturbative for typ- 
ical values of the lattice spacing. They can therefore be 
removed systematically using the method of Symanzik. 
This gives lattice operators with better taste symmetry 
and smaller scaling violations. We do this by smear- 
ing the gauge field to apply a form factor to the quark- 
gluon vertex and suppress the coupling between quarks 
and high-momentum gluons |39l |4Q, 41]. This amounts 
to replacing the gauge link U in Eq. (|17() by a "fat" link: 
a weighted sum of staples as shown graphically in Fig.^ 
The choice of weight factors, determines the degree of 
improvement. In this paper we compare three different, 
improved operators known as AsQTAD, FatTxAsq and 
Hyp. 

The first two are based on the Fat 7 operator which 
includes staples of length up to 7a. We can remove all 
0(0^) discretization errors (at tree level) by adding two 
more terms to the Fat7 prescription, including a Naik 
three-link term and a five- link Lepage staple [s^ . 
This operator is called AsQ, whereas Asqtad is the tad- 
pole improved version which aims to also reduce the ra- 
diative corrections to the scaling. 

The improvement procedure may be iterated; for our 
most highly improved operator the original links are first 
Fat7 smeared and projected back onto the SU (3) gauge 
group. These fattened links are then used to construct 
the AsQ operator. We call this operator Fat7xAsq. We 
stress that this is the same operator that was used in . 
Regrettably, there it was misnamed HiSQ, a name already 
used to describe a slightly different Dirac operator. 




e3' (Naik) 



FIG. 1; Paths entering the improved staggered operators. 



The final improved staggered Dirac operator is Hyp 
|43| . Its construction is motivated by perfect action ideas, 
and involves three levels of (restricted) APE smearing 
with projection onto SU (3) at each level. The restrictions 
are such that each fat link includes contributions only 
from thin links belonging to hypercubes attached to the 
original link. 

The spectrum of all improved Dirac operators is of the 
form in Eq. 119(1 . Improved operators, however, show 
very small taste-changing effects both in the hadronic 
splittings [23, 01 and, as we shall see, the eigenvalue 
spectrum. 



III. DETAILS OF THE SIMULATION 

In this section we provide details of the gauge back- 
grounds used in this study, and the methods used to de- 
termine the eigenvalues and the topological charge. We 
also introduce notation for the eigenvalues to be used in 
the following sections. 



A. The gauge action 

We use a quenched, SU (3) gluonic action that is both 
tree level Symanzik and tadpole improved pslli^li^liM 

(23) 

where P, R are 1x1 and 2x1 Wilson loops respectively. 
The tadpole improvement coefficient uq is defined as the 
fourth root of the mean plaquette. 

We have generated five ensembles of around 1000 con- 
figurations, with parameters for these ensembles as in 
Table ^ (taken from Refs. [sol Isij ) . The physical spatial 
size of the lattice is aL, and the volume a'^V = a^L^T. 
We have also used a few configurations from a Wilson 
gauge ensemble with lattice spacing a « 0.1 fm and vol- 
ume 16^ X 32. 

Three ensembles have been tuned to study the effect 
of varying the lattice spacing at fixed physical volume 
(aL « 1.5 fm). Another three vary the volume at fixed 
a = 0.093 fm. 

The two coarser lattice spacings are representative of 
current, large scale unquenched simulations using im- 
proved staggered quarks, with the smallest a being in- 
dicative of future calculations. Our gauge action differs 
from the one used by the MILC Collaboration js^ only in 
small one loop radiative corrections which do not affect 
the relevance of the results in this paper. 

The main objective of this paper is to compare our 
results with the expectations from the continuum limit 
and with the predictions from RMT. For this purpose any 
systematic uncertainty in the determination of the lattice 
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TABLE I: Simulation parameters: the gauge coupling /3, lat- 
tice spacing and volume and the spatial lattice extent. 



J3 g/fm V = L^T aL/im 

4.6 0.125 12* 1.50 

4.8 0.093 12* 1.12 

4.8 0.093 16* 1.49 

4.8 0.093 20* 1.86 

5.0 0.077 20* 1.54 



spacings (for example due to finite-volume effects) is not 
important. It would only become important for obtain- 
ing a precise determination of some physical quantity. 
Therefore in general we ignore this source of errors in 
our measurements, using only the statistical errors based 
on a jack-knife analysis of the data. 



B. Measuring the topological charge 

We measure the topological charge Q using a number 
of gluonic methods to look for consistency. The topo- 
logical charge of a gluonic configuration is not uniquely 
defined at finite lattice spacing, and different methods 
will sometimes disagree. Over an ensemble, the number 
of configurations that disagree on the value of Q vanishes 
as a? in the continuum limit. 

We first cool the gauge fields and then apply a lattice 
topological charge operator. The cooling is (separately) 
done with two different gauge actions (the Wilson and the 
5Li 01 )• We measure the topological charge with 
the hig hly accurate "5 Loop improved" (5Li) operator 
[s^. l54j . We then check that both charges are very close 
to integers and stable under cooling. If the topological 
charges (rounded to the nearest integer) defined by this 
procedure are the same for both actions, then we assign 
the configuration to a definite topological charge sector. 

Does this consistency criterion bias our fixed-Q 
subensemble averages? No; at a = 0.093 fm ensembles 
fewer than 10% of the configurations have an ambiguous 
topological charge. For the finest ensemble, a = 0.077 fm, 
this number goes down to around 2%. 

We have repeated parts of the analysis by ignoring the 
consistency requirement and using only one of the cooling 
actions in the definition of Q, assigning every configura- 
tion to a topological charge sector. The results are in- 
distinguishable from those presented here, and therefore 
this robustness criterion neither introduces a systematic 
uncertainty in our results, nor changes the conclusions. 

We stress that we only use cooling to determine the 
topological charge. All the Dirac spectrum measure- 
ments are done on the original thermalized ( "hot" ) con- 
figurations. 



C. Determining the eigenvalues 

To calculate the low-lying spectrum of the Dirac opera- 
tors we use an implementation of the CuUum-Willoughby 
Lanczos-based algorithm [s^ ■ Due to the structure of the 
staggered Dirac operator the calculation of the spectrum 
can be done on an hermitian, positive semi-definite oper- 
ator, defined only on half of the lattice (either the odd or 
the even sublattice). The algorithm can be extended to 
the calculation of the corresponding eigenvectors, which 
are needed to obtain the chirality. 

On the ensembles used in this study we found no signif- 
icant numerical difficulties with this algorithm. In prac- 
tice, the absence of exact zero modes simplifies the cal- 
culation from a numerical point of view, as the linear 
operator we are dealing with is then positive definite. 

Due to the exact ± symmetry, we only need to consider 
half of the spectrum, say {lAg, As > 0} (which we call the 
positive sector) . The other half is a mirror image which 
we call the negative sector. From now on we dispense 
with the imaginary unit i, and discuss and show on the 
plots the imaginary part of the positive half spectrum 
{As, As > 0}. 

For the analysis we divide the (positive) eigenmodes 
on each configuration into two sets based on the gluonic 
topological charge. Inspired by our continuum limit ex- 
pectations, we expect the 4|(5| near-zero modes to divide 
into 2|Q| modes each side of zero. We call the arithmetic 
mean of the positive near-zero modes Ap. 

The remaining positive modes are divided into quar- 
tets. We call the means of successive sets Ai^2,...- We de- 
fine the intra-quartet splitting of the s'^ set, 5 Kg, as the 
difference between the largest and smallest eigenvalues in 
the quartet. We also define the inter-quartet gap between 
the s*** and {s + 1)*^ quartets, AAs, as the difference be- 
tween the smallest eigenvalue of the higher quartet and 
the largest eigenvalue of the lower quartet. 

We use the notation ( • )q to denote the expectation 
value over sub-ensembles with fixed gluonic topological 
charge itQ. 

Unless otherwise specified, the eigenvalues are ex- 
pressed in physical units of MeV, with the scale set using 
the quoted lattice spacing. 

IV. SPECTRUM AND INDEX THEOREM 

In this section we study the Dirac spectra, and show 
how the continuum results emerge for the improved op- 
erators. 



A. The Index Theorem 

We begin our analysis by qualitatively comparing the 
low-lying modes of various staggered quark operators 
on typical gauge backgrounds selected from those with 
Q = 2. Working at approximately fixed lattice volume 
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FIG. 2: The positive half of several typical low-lying eigenmode spectra for configurations of Q = 2, in lattice units. The 
X-axis is eigenvalue number. The Hyp operator gives results very similar to Fat7xAsq, and they are not shown for clarity. 



(aL « 1.50 fm), we compare three different gauge field 
actions: unimproved Wilson at a = 0.1 fm; Symanzik im- 
proved, Eq. (|23|l . at a comparable a — 0.093 fm, and at 
the finer a — 0.077 fm. In Fig.|2|we show both the value 
and the chirality of the first sixteen positive eigenvalues. 

Near the continuum limit we expect to see first 
2\Q\ — 4 near-zero modes with their chirality renor- 
malised slightly away from unity. The remaining modes 
should have chirality near zero and divide into almost 
degenerate quartets with intra-quartet splitting much 
smaller than the inter-quartet gap. 

For the improved gauge action on the finer lattice the 
spectrum looks quite continuum-like for all Dirac opera- 
tors and we see a clear Index Theorem. For the improved 
operators we also see a very clear quadruple degeneracy 
in the non-zero modes. The renormalisation of the chi- 
rality away from 1 is small for the improved Dirac op- 
erators (around Z = 1.2), becoming as large as Z ^ 4 
for the One-link fermions. This result is fairly typical 



of operator renormalisation for staggered quarks, with 
the improved actions sho wing significantly better renor- 
malisation factors |5^ [s^, Eg. We will discuss later the 
renormalisation of the chiral condensate and see the same 
effect. 

On the coarser improved gauge background, we still 
have a good approximation to the Index Theorem using 
the FatTxAsq and Hyp operators: there is a separation 
between near-zero modes of high chirality and the rest 
of the spectrum; the number of near-zero modes is 2\Q\; 
and all of them have very similar chiralities. The rest 
of the spectrum have almost zero chirality and cluster 
into quartets. These properties can be seen to a lesser 
extent for the AsQTAD operator, but are absent in the 
One-link case. 

Nothing clear can be seen in the unimproved, Wilson 
gauge case. Improving the Dirac operator makes a differ- 
ence, increasing the chirality of the low modes. Nonethe- 
less the chiralities do not show a clear Index Theorem, 
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FIG. 4: A scatter plot for Fat7xAsq, with absolute value of 
the chirality on the y axis and eigenvalue As on the x axis, in 
lattice units. The lowest 50 eigenvalues for around 100 con- 
figurations from our 20*, a = 0.077 fm ensemble are plotted. 
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FIG. 3: A scatter plot for different staggered quark formu- 
lations, with absolute value of the chirality on the y axis and 
eigenvalue As on the x axis, in lattice units. The lowest 50 
eigenvalues for 147 configurations from our 16", a=0.093 fm 
ensemble are plotted. 



and the eigenvalues do not gather into quartets. 

To give an idea of the behaviour over a large ensemble. 
Fig. U] shows a scatter plot of the absolute value of the 
chirality versus the (positive) eigenvalues. Results for dif- 
ferent operators are plotted, using the same a = 0.093 fm 
ensemble of gauge fields. We can see the formation of a 
gap between modes of small and large chirality as we 
improve the staggered operators, as well as the overall 
increase in chirality of the near-zero modes. As seen in 
Fig. ^ for the most improved operator on the finest en- 
semble, the separation gets sharper as we approach the 
continuum limit. 

The separation is not strict even for the most improved 
action, with a few configurations showing intermediate 
values of the chirality. At least some of these cases are as- 
sociated with configurations where the 0{a^) ambiguity 
in the topological charge is large and/or where the total 
number of near-zero modes is greater than the minimum 
prescribed by the Index Theorem, i.e. n+ +n_ > A^tlQ]. 

As the near-zero modes for the more improved actions 
clearly separate and have well defined chirality, we may 
define a fermionic index Q = (n+ — n-)/Nt on each con- 
figuration. Q is then strongly correlated with Q, and we 
have a very good approximation to the Index Theorem. 
We expect the number of configurations where Q ^ Q to 
vary as 0{a'^). Although we don't have enough statistics 
to test this scaling quantitatively, we can get some esti- 



mates from our results. For our a = 0.093 fm ensemble, 
the disagreement between Q and Q is around 10%, and 
it goes down to about 2% on the finest a = 0.077 fm 
ensemble. This ambiguity is therefore of the same order 
as the one intrinsic in the gluonic definition, as discussed 
in Sec, nil I3l We stress here that Q played no part in the 
assignment of configurations into sectors of fixed topo- 
logical charge. 

We may worry that when Q ^ Q the division of the 
eigenvalues into near-zero and non-zero modes, and the 
division of the latter into quartets, will be skewed. For 
instance, on such configurations the intra-quartet split- 
ting will be as large as the inter-quartet gap. The exam- 
ination of the histograms of SAs and AAg show that in 
all cases any slight bias on the mean is swamped by the 
statistical spread. 



B. The Banks— Casher relation 

We begin our quantitative analysis with the Banks- 
Casher relation (|lll) . which links the density of non-zero 
eigenvalues at the origin with the chiral condensate, and 
thus provides a model independent determination of S. 
For Nt degenerate fiavours, the Banks-Casher relation 
would be modified to 



S = np{0)/Nt 



(24) 



This is the relevant form for staggered fermions, with 
Nt — 4. Note that the appropriate ijjip operator for the 
chiral condensate is the taste-singlet, but this is the local 
operator and so there are no other modifications to the 
Banks-Casher relation for staggered quarks. 

In Fig. |31 we plot a histogram of the spectral density 
averaged over all configurations in the ensemble. The 
errors come from a jack-knife analysis of the data. 

The spectrum divides clearly into two separate sec- 
tions: a sharply peaked spike near A = 0, coming from 
the near-zero modes, and a broader distribution at larger 
eigenvalue from the remaining modes. This separation is 
much sharper for the improved actions. The spectral den- 
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TABLE II: Estimates of the cube root of the bare chiral 
condensate, E^''^, (in units of MeV) from: the Banks-Casher 
relation (upper section) and comparisons of eigenvalue distri- 
butions with random matrix theory (lower section). 

Dirac operator 



a/fm 


aL /fm 


One-link 


ASQTAD 


Hyp 


Fat7xAsq 


0.125 


1.50 


348 (10) 


274 (10) 




256 (10) 


0.093 


1.12 


330 (30) 






260 (10) 


0.093 


1.49 


349 (10) 


275 (10) 


279 (10) 


260 (10) 


0.093 


1.86 


360 (10) 






255 (10) 


0.077 


1.54 


322 (10) 


253 (10) 




238 (10) 


0.125 


1.50 


356 (8) 


283 (6) 




260 (6) 


0.093 


1.12 


383 (43) 






289 (33) 


0.093 


1.49 


361 (12) 


293 (5) 


279 (4) 


281 (2) 


0.093 


1.86 


358 (6) 






286 (4) 


0.077 


1.54 


321 (10) 


266 (10) 




256 (5) 



sities corresponding to the improved actions also show an 
excellent scaling behaviour with the lattice spacing. 

For larger A, a A"^ behaviour is expected to set in. Hav- 
ing only calculated the first 100 eigenvalues, however, we 
instead find the spectral density cut olT at larger A. It 
would be necessary to include many more modes to see 
this behaviour. 

Using Eq. H24() . we can extract the single taste chi- 
ral condensate by extrapolating the density of non-zero 
modes to zero, p(A) p(0). The chiral condensate is 
the order parameter for chiral symmetry breaking, and 
we see clear evidence from Fig. [S] that all our lattices are 
large enough that E ^ 0. 

To estimate p(0) we carry out uncorrelated fits to p(A), 
excluding the near-zero mode peak from the fit range. We 
use a polynomial fit function 



p{\) = po + Vi + \^P2 + >?Pz 



(25) 



We carry out a number of fits to each data set, constrain- 
ing different combinations of pi to zero and varying the 
fit range. The intercept po is then used to obtain the chi- 
ral condensate values given in the upper half of Table ITU 
The quoted errors cover the fit function and range uncer- 
tainties, which dominate the statistical and fitting errors. 

Note that the One-link density at the origin p(0) is 
much larger than the one for the improved case. This 
is clearly noticeable when inverting the operator, with 
the improved actions needing fewer iterations than the 
One-link action (also seen in [ilj). 

The larger p(0) gives a much higher value for S for 
the One-link action, but it is important to note that 
this is the bare lattice chiral condensate and does not 
have to agree between actions in the same way that the 
bare lattice quark mass does not. What should agree 
is the chiral condensate defined in some physical way, 
for example in the MS scheme at a fixed scale. The 
appropriate renormalisation factor for tpij) is the inverse 
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FIG. 5: The spectral density for different Dirac operators at 
fixed lattice volume aL = 1.5 fm. Only the contribution of 
the first 100 eigenvalues have been calculated. 



of that for the quark mass and this has been calculated 
to 0{as) for the staggered quark actions used here and 
the improved gluon action (at 0{as) there is no effect 
from one-loop improvements to the gluon action or the 
presence of sea quarks). The quark mass renormalisation 
was quoted in 59] as: 

m^{p) = 

^ {l + av{q*) Z(^\ap,{am)o) + 0{a^)) , 

Z^^\ap,amo) = (^b{amn) - - ^\n{ap)^ . (26) 

where 6(ar7io) is the finite piece of the lattice mass renor- 
malisation and has small anio dependence. The comes 
from the continuum mass renormalisation. With tadpole- 
improvement b{0) takes the value 0.543 for the Asqtad 
quark action and 3.6 for the One-link action (denoted 
naive staggered quarks in jH^)- For FatTxAsq b{0) is 
0.375. For the chiral condensate renormalisation, Zs, we 
need to take the inverse of the renormalisation above. 
We take to be 2 GeV, as a standard reference point, 
and take the scale of to be 2/a, as determined for the 
Asqtad action at comparable values of a in The 
values of as in the V scheme in the quenched approx- 
imation are taken from 60]. It is then clear that Zs 
for the Asqtad and FatTxAsq actions is very close to 
1 for all lattice spacing values, but significantly smaller 
than 1 for the One-link action. The renormalisation 
brings the value of the physical chiral condensate for the 
One-link action from the Banks-Casher relation down 
to 290 MeV, much closer to that from the improved ac- 
tions. The renormalisation constant is so far from 1 in 
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this case, however, that significant higher order correc- 
tions are to be expected. For the improved actions, how- 
ever, the renormaUsation is very small indeed and num- 
bers change by less than 2% from the bare results of 
table m The estimated uncertainty from higher order 
perturbative corrections is 4% i.e. 1 x a^. From table Hll 
we also have a statistical error of 4%. In addition there 
are statistical and systematic errors from the determi- 
nation of the lattice spacing that are relevant here. All 
of these errors are small compared to the overall 10-20% 
error from the inability to determine the lattice spacing 
uniquely in the quenched approximation. This means 
that any physical value for the chiral condensate that 
we quote will have a significant error from this source 
alone. Our results for the FatTxAsq action give a re- 

suh for {-{^^PY'^Y''^ (2GeV) of 260(10)(10)(30) MeV, 
where the first error is statistical, the second perturba- 
tive and the third a quenching error. This is in agreement 
with other results in the literature in the quenched ap- 
proximation (for example, HJ). An improvement on this 
value requires unquenched configurations, but our result 
shows that an accurate calculation is certainly possible 
with improved staggered quarks. 



C. Scaling of the near-zero modes 
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FIG. 6: The scaling of the near-zero modes with lattice spac- 
ing at fixed lattice volume for different topological charge sec- 
tors. The near-zero modes are averaged. Also shown are lin- 
ear plus quadratic fits through the origin for Asqtad and 
Fat7xAsq . 



The near-zero eigenvalues for a given Dirac operator 
vary both with the lattice spacing and volume. 

In Fig. we show the variation of the mean near-zero 
modes (Aq) with lattice spacing at fixed physical vol- 
ume [ah w 1.5 fm). The panels show the results for 
sectors of topological charge = 1 and 2. For a suf- 
ficiently small lattice spacing, we expect the eigenvalues 
to go to zero as {a}) plus 0{a^^ corrections. The im- 
proved staggered results show this trend across the full 
range of lattice spacings studied, with small 0{a^') cor- 
rections, slightly positive for FatTxAsq but negative for 
Asqtad. By contrast, the One-link near-zero modes 
show very strong deviations from the leading order scal- 
ing behaviour. 

In almost all quantities that we measure the Hyp oper- 
ator gives results that agree with FatTxAsq within one 
standard error. Unless there is a significant difference 
we do not discuss the Hyp data and for clarity they are 
not shown on the figures, but we have checked that all 
statements regarding FatTxAsq apply equally to Hyp. 

The FatTxAsq near-zero modes are numerically 
smaller than those of the One-link action for a < 
0.13 fm. The Asqtad near-zero modes are smaller than 
the One-link ones for a < 0.08 fm. Even when the 
modes are similarly sized, however, we have seen earlier 
that there is a qualitative difference in the chirality as 
we improve the action. This is a clear signal that the 
operator improvement is working. 

The variation of the near-zero eigenvalues with vol- 
ume at fixed lattice spacing is shown in Fig. [7| (Aq) is 
expected to vary as the inverse square root of the volume. 



This can be justified in a semiclassical picture, where the 
topological charge density is dominated by a dilute gas 
of instantons. The finite volume corrections to the zero 
modes come from the lack of large instantons, excluded 
from the lattice by their cores sizes r > aL. Such cor- 
rections should vary as or l/VaH^ |62l |. 

We see from Fig.|2|that the data at fixed lattice spacing 
follows this approximate volume scaling quite closely. We 
might ask whether (Aq) really goes to zero in the large 
volume limit at finite lattice spacing. Fits to our data 
allowing a non-zero intercept are not stable, but a large 
reduction in statistical errors, as well as simulations in 
larger volumes, would be needed to answer this question 
definitively. 



D. Scaling of the low- lying non-zero modes 

We turn now to the low-lying non-zero modes. From 
the requirement of a physical spectral density in Eq. I|10|) 
and the Banks-Casher relation, Eq. (|ll|l . it is clear that 
the size of these is governed, up to a constant, by the 
product of the bare chiral condensate and the volume 
TiV. Specifically, the number of modes in a small inter- 
val (0, S) should scale as ^ 6p{0)V ^ SJ^V, and therefore 
we expect each mode to scale as (Y,V)~^. The variation 
of these modes with changing lattice spacing and volume 
is therefore fixed by the response of the bare chiral con- 
densate. 

In Fig. 13 we show the variation of the first two non- 
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FIG. 7: The averaged near-zero modes as a function of lat- 
tice volume at a = 0.093 fm for different topological charge 
sectors. Also shown are linear plus quadratic fits through the 
origin. 
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FIG. 8; The scaling of the non-zero modes with lattice spac- 
ing at fixed lattice volume for different topological charge sec- 
tors. In each case the quartets are averaged. 



zero mode quartets (Ai_2) with lattice spacing at fixed 
physical volume (aL k, 1.5 fm) for several Q. For all the 
improved actions the changes are very small, indicating 
only small discretization errors, as for the chiral conden- 
sate. For the near-zero modes, the deviations from the 
leading order scaling behaviour seem to be of opposite 
sign for AsQTAD and FatTxAsq. The One-link data 
show much larger deviations, of the same sign as for the 
AsQTAD data. 

The different actions appear to lack a common contin- 
uum limit. This is not surprising: as discussed above. 



the plot is inversely related to the density and therefore 
the bare chiral condensate. As explained earlier, this re- 
ceives a large renormalisation for the One-link action 
compared to the improved actions that means that the 
bare results are not comparable. 

From the arguments above, we expect (As) to vary as 
the inverse volume at fixed lattice spacing. We find this 
to be so for both One-link and improved actions. Fig.|2| 
shows the results for s = 1, 2. 

Given the different volume scaling of the two sets of 
eigenmodes, we cannot really talk about a gap between 
near-zero and non-zero modes: we can always choose a 
volume large enough that the "non-zero" modes (sup- 
pressed as 1/1^) are numerically smaller than the "near- 
zero" modes (which fall only as Xj^/V). The modes are 
still qualitatively different however, and can be separated 
by their chirality. 

Taste breaking interactions disappear in the continuum 
limit and the eigenmodes should recover a 4-fold degen- 
eracy. We plot the intra-quartet splitting of the first 
quartet, ((5Ai) for different sectors of topological charge 
in Fig. [TUl fhigher quartets have smaller splittings). The 
lattice volume is fixed at ah « 1.5 fm. The splitting 
falls with lattice spacing, as expected, and in a way con- 
sistent with the 0{a^) dependence of the taste breaking 
interactions. The splitting depends only weakly on the 
topological charge. 

This is the only quantity we have studied where there 
is a difference of more than one standard error between 
Hyp and FatTxAsq results; the FatTxAsq results are 
marginally closer to our continuum expectations. (This 
was also true in Fig.|Bl but the difference was not statis- 
tically significant). 

There are strong similarities between this plot and 
Fig. El the improved fermion actions all show only very 
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FIG. 12: The ratio of the intra-quartet splitting of the first 
non-zero quartet to the inter-quartet gap between first and 
second quartets as a function of lattice spacing at fixed lattice 
size aL ~ 1.5 fm. The lines are linear plus quadratic fits for 
Asqtad and Fat7xAsq . 
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FIG. 11: The scaling of the first non-zero intra-quartet split- 
ting with volume at fixed lattice spacing in different topolog- 
ical charge sectors. Higher quartets have smaller splittings. 
See comments in Fig. |S| 



To satisfy ourselves that there is a clear restoration of 
the taste symmetry, we need the intra-quartet splitting 
((5As) to be much less than the gap between neighbouring 
quartets (AA^). In Fig.^]we show the ratio of these as a 
function of lattice spacing for s = 1 (where it is largest) . 

Again, the ratios show little dependence on the topo- 
logical charge. In all cases the improved action results 
are better than the One-link. This is increasingly true 
at small lattice spacings, as the improved actions show a 
ratio falling clearly to zero as a? . 

We can understand these results quite easily. As the 
quartets become degenerate, (AAg) will tend towards the 
difference between the quartet means. Fig. |S1 shows this 
is controlled by the chiral condensate and has little lattice 
spacing dependence for the improved actions. The intra- 
quartet splitting falls as a? and so, therefore, will the 
ratio. 



weak deviations from the naive scaling behaviour; the 
O(a^) corrections are positive for FatTxAsq and neg- 
ative for Asqtad; the One-link results however show 
very large deviations from the expected behaviour. Nu- 
merically, the FatTxAsq splittings are smaller than the 
One-link for a < 0.12 fm; the Asqtad numbers are less 
for a < 0.08 fm. 

In Fig. 111! we show the variation of the splitting with 
volume at fixed lattice spacing. Once more the splitting 
is independent of the topological charge. There is a strik- 
ing difference in the volume dependence: the One-link 
splittings vary markedly whilst the improved action re- 
sults are nearly insensitive. 



V. RMT PREDICTIONS 

As discussed previously, the chiral symmetries of stag- 
gered quarks are more complicated than for continuum 
QCD. At finite lattice spacing the = 16 pions split 
into 5 multiplets containing (1, 4, 6, 4, 1) states, and only 
one of these states becomes massless in the chiral limit. 
This is known as the "Goldstone pion", with mass Mg- 
The fifteen remaining states have masses -Mng arising 
from the taste breaking interactions. The (chirally ex- 
trapolated) masses are 0{a^) and therefore zero only in 
the continuum limit. 

There are potentially two universal regions for different 
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FIG. 13: The ratios of expectation values of small eigenvalues 
for the aL = 1.5 fm, a = 0.093 fm ensemble for different Dirac 
operators, compared with the predictions based on a universal 
distribution (horizontal lines) for IQI < 2. 
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FIG. 14: The ratios of expectation values of small eigenval- 
ues for the aL = 1.5 fm ensembles for the Fat7xAsq operator 
compared with the predictions based on a universal distribu- 
tion (horizontal lines) for |Q| < 2 . 

parameters of the system: 

e-regime: (Aqcd)"^ < i < (Mng)"^ , 
e'-regime: (A/ng)"^ < i < {Mq)-^ . 

The first corresponds to Eq. \i2\ . with the same chiral 
symmetries that we expect in the continuum. The low- 
lying non-zero modes have a near A'^^-fold degeneracy, and 
follow the predictions of Eqs. H14I15|I . These are exactly 
the same distributions as for continuum QCD (and chiral 
lattice fermions |H1 El El ElEI ) . 

At finite lattice spacing there is a second, e'-regime, 
where the finite volume partition function describes only 
the static mode of the single Goldstone pion. There is not 



even approximate restoration of the continuum symme- 
tries, and the associated RMT has only U (l)(g)J7(l) chiral 
symmetry. The universal predictions for the eigenvalues 
are therefore strikingly different from those for contin- 
uum QCD. In particular, the predictions are the same 
for all sectors of topological charge [63| . Presumably it is 
this regime that was studied in llllillllllllll. With 
the coarse lattices and unimproved gauge ensembles used 
in these studies, they observed no sensitivity of the eigen- 
value spectra to Q, leading to the incorrect folklore that 
staggered quarks are "blind to the topology" . 

To see the continuum chiral symmetries, there must be 
a sufficiently large mass gap between the heaviest pion 
and the lightest of the other hadrons. This requires the 
use of improved gauge and fermion actions and a suffi- 
ciently small lattice spacing. In particular if we increase 
the size of the lattice L, the lattice spacing must be re- 
duced accordingly in order to satisfy L <^ (Mng)~^. 
Otherwise we would start to see the effects of the mass 
of the non- Goldstone pions. At coarse lattice spacings, 
with unimproved fermions, A/ng ~ Aqcd and there will 
be no e-regime. 

If the lattice spacing is keep fixed and we increase the 
volume, we expect there will be a non-universal crossover 
until the lattice reaches the e'-regime when universal be- 
haviour will resume. Of course, it would be interesting to 
simulate a range of lattice sizes and witness the crossover 
between the two universal regions. The range of lattice 
spacings is, however, well beyond the resources of this 
study. We instead concentrate on the continuum-like re- 
gion. 



A. Comparison 

We begin by comparing the ratios of the mean non- 
zero eigenvalues {As)q / {At)Q (which we denote "s/t") 
to the ratios from RMT. Factors of the chiral condensate 
scale away, so this provides a parameter-free test of the 
e-regime predictions. 

We show the eigenvalue ratios in Figs. I13I15I The ra- 
tios differ markedly between different topological charge 
sectors, showing clearly that the staggered fermions are 
in no way "topology blind". Also shown are the ra- 
tios predicted by RMT. Fig. 1131 shows that all the stag- 
gered actions agree quite well with the expected ratios at 
a = 0.093 fm. We therefore have good evidence for the 
existence of a continuum e-regime. It also implies that 
the continuum chiral symmetries are present, with the 
correct pattern of spontaneous breaking. 

There are some deviations from the universal predic- 
tions. They are not discretization effects: Fig. 1141 shows 
the FatTxAsq ratios for a range of different lattice spac- 
ings at fixed volume aL « 1.5 fm and discrepancies exist 
even on the finest lattice. Increasing the system size from 
1.5 to 1.8 fm in Fig. E| does reduce the size of the dis- 
agreement. Indeed, on the larger volume the ratios agree 
with RMT at least as well as for chiral fermions on com- 
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FIG. 15: The ratios of expectation values of small eigenvalues 
of the Fat7xAsq operator for two a = 0.093 fm ensembles of 
different volumes compared with the predictions based on a 
universal distribution (horizontal lines) for \Q\ < 2. 

parable lattices [68l |. 

Fig. El appears to suggest that the One-link Dirac 
operator fits the universal predictions as well as the im- 
proved operators. This appears so even on the coarsest 
lattices despite the large pion splitting and lack of obvi- 
ous spectral quartets in the One-link spectrum. 

In this respect the "ratios plots" can be misleading. 
Firstly they only compare the mean eigenvalues with 
RMT and provide no information as to whether the shape 
of the individual distributions match. Also, when com- 
paring with RMT we must rescale all the eigenvalues us- 
ing the chiral condensate. Agreement with the ratio plot 
does not imply that eigenvalues in different topological 
charge sectors are scaled by the same factor. 

To address the first point we separately fit the individ- 
ual spectral densities to the predictions from RMT. These 
one parameter fits based on Eq. 1)13(1 yield a prediction 
for the chiral condensate. To avoid the ambiguities of 
histogram bin-size we use the cumulative eigenvalue dis- 
tributions, 

p{X) = I dX' p(A') . (27) 
Jo 

We show examples of the fits in Figs. 1161 and 1171 

For all operators the fits are worst for low s and for 
small \Q\. We can qualitatively understand this in terms 
of eigenvalue repulsion. The taste changing interactions 
have split the near-zero modes, which has a knock-on ef- 
fect on the rest of the spectrum as the eigenvalues avoid 
level crossings. The lower quartets are closest to the near- 
zero modes and are thus most strongly affected. As the 
quartets get numerically larger as we increase the topo- 
logical charge, so the repulsion should be less. 

It is clear that the improved Dirac operators match 
the predictions of RMT very closely over the full range 



of lattice spacings. By contrast, even at a = 0.077 fm the 
One-link Dirac operator still shows significant discrep- 
ancies. 

Separately fitting the {Q — 0,k = 1..4) and {\Q\ ~ 
0..3, fc = 1) eigenvalue distributions to the RMT predic- 
tions gives 7 values for the bare chiral condensate for each 
Dirac operator on each ensemble. We find these numbers 
to be broadly consistent. The spread in the numbers is 
greatest for the One-link action and on the smallest vol- 
umes. This is not surprising, when such eigenvalue dis- 
tributions showed the greatest deviation from the RMT 
predictions 

The lower half of Table HTl gives estimates for the chiral 
condensate obtained from the median of the seven indi- 
vidual numbers. The quoted errors are half the range, 
which forms the dominant uncertainty. 

We can compare these with the model-independent es- 
timates of S obtained using the Banks-Casher relation. 
There is very good agreement between the numbers. The 
dominant source of error is different in each case: for 
the Banks-Casher relation it comes from the choice of 
fit range and function used to extrapolate the spectral 
density to zero eigenvalue. For the RMT comparisons, it 
comes from the variation in fit parameters for different 
eigenvalues, and is therefore related to exactly how well 
the effective partition function describes the lattice the- 
ory. The good agreement of the numbers suggests both 
systematic biases are under control. See section ITV Bl for 
a discussion of how the bare chiral condensate must be 
renormalised to extract a physical result. 



VI. CONCLUSIONS AND OUTLOOK 

It is widely believed that any "problems" with a given 
lattice Dirac operator will show themselves most clearly 
in the topologically sensitive, low-lying eigenmodes. In 
this study we have studied these modes for the staggered 
lattice Dirac operator, using a range of improved and 
unimproved versions. The improved operators pass all 
the tests. In particular, and contrary to previous ac- 
cepted wisdom, such fermions do respond exactly as ex- 
pected to the gluonic topological charge, and in a way 
identical to continuum QCD and other lattice formula- 
tions of the Dirac operator. 

We have seen that for improved operators the eigen- 
value spectrum divides cleanly into near-zero and non- 
zero modes. The near-zero modes are characterised by a 
uniformly high chirality, with their relative number fixed 
by the Atiyah-Singer Index Theorem. Indeed they could 
be used to define the topological charge Q of a configu- 
ration. Their eigenvalues scale with lattice spacing and 
volume as expected. 

The non-zero modes, by contrast, have chirality that 
is near zero. They divide into near degenerate quartets, 
with the splitting reducing to zero as in common with 
the taste breaking interactions. The quartet means are 
controlled by the chiral condensate and scale with lattice 
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FIG. 16: Comparison of the unfolded eigenvalue distribution with RMT for the lowest eigenvalues in the Q = sector. The 
lattice size is kept at aL ~ 1.5 fm 



spacing and volume accordingly. In addition, the low- 
lying non-zero modes follow closely the universal distri- 
butions predicted by random matrix theory. 

Since the near-zero and non-zero modes scale differ- 
ently with the volume, it is clear that the "gap" between 



the highest near-zero and the lowest non-zero modes is a 
volume dependent quantity. At finite lattice spacing we 
can always go to a volume large enough that the two sets 
of eigenvalues have comparable magnitude. This is not 
really an issue because the gap is not a physically observ- 
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FIG. 17: Comparison of the unfolded eigenvalue distribution with RMT for the lowest eigenvalue for different topological 
charges. The lattice size is kept at aL ~ 1.5 fm 



able quantity: what matters is the size of the near-zero 
modes compared to the light quark masses and whether 
there is a well-defined Index Theorem. As we have seen 
the two sets of modes are distinguishable by their chiral- 
ity, even in the large volume limit. 



A quantitative study of axial anomaly physics will re- 
quire the near-zero modes to be smaller than the light 
quark masses, i.e. A < 0(5MeV), or, at least, smaller 
than the sea and valence quark masses used in the simu- 
lation. Then the eigenvalue spectrum is cut off at the low 
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end by the quark mass and the fact that the near-zero 
modes are not actually zero is irrelevant. 

To understand what this means for our quenched en- 
sembles, it is necessary to define the size of the typical 
near-zero mode. The mean zero mode (Ao)q increases 
as we increase Q, due to eigenvalue repulsion. For our 
f3 = 4.8, L = 20 lattice we find that the mean zero- 
mode is fitted quite well by (Ao)q « (3 -I- Q) MeV for 
Q = 1..3. If we define the typical topological charge Q 
to be the RMS value from the topological susceptibil- 
ity, Q ~ VxV — 2.6, then the typical near-zero mode 
is Ao = (3 2.6) w 6 MeV. This suggests that we are 
already able to study the quenched axial anomaly. 

We can attempt to extrapolate our quenched results 
to ask whether this condition is true on existing con- 
figurations from the MILC collaboration This re- 
quires various unsupported, if reasonable, assumptions, 
not least that sea quark effects do not strongly affect the 
size of individual near-zero modes and that (Ao)q con- 
tinues to be linear in Q. 

Typical dynamical simulation volumes are larger than 
we have used in this quenched study, to avoid finite vol- 
ume effects in hadron spectroscopy. As well as an over- 
all Vv suppression in the size of the near-zero modes, 
this will also affect the typical value Q and therefore the 
amount of eigenvalue repulsion. There are two compet- 
ing effects: the sea quark induced screening of topological 
charge reduces x, whilst the larger volume increases the 
width of the histogram of Q. 

Consider the a — 0.09 fm MILC ensemble with sea 
quark masses nis/mu d = 5 and volume V = 28'^ x 96. 
Using the data in z^, we obtain Q = 7.3. As- 
suming the quenched relation above holds, we have 
Ao = (3 + 7.3)/v/284/204 - 6 MeV (to be conservative 
we assume that the volume scaling of the near-zero modes 
goes with the smallest length of the lattice, therefore 28^ 
instead of 28^x96). This is very encouraging. The light- 
est sea quark mass presently used in these dynamical sim- 
ulations is mo = 16 MeV, which is several times larger 



than this projected near-zero mode. 

Even given the caveats surrounding this rough esti- 
mate, it does seem likely that a systematic study of the 
QCD axial anomaly is possible on these configurations 
and various aspects of this work are now in progress. 

An on-going debate surrounds the theoretical justifi- 
cation of the methods use to simulate Nf < Nt fiavours 
of sea quarks rZlj. It^ ItJ . Recent work suggests there 
is no ambiguity [79, UM- Our results (and related 
work elsewhere |3.'tI 123) indicates that with improved 
staggered quarks the spectrum matches very closely the 
continuum expectations with clear restoration of taste 
symmetry. It therefore paves the way for the full analy- 
sis of the staggered taste basis that is necessary to finally 
resolve this issue. 

To clarify the role of topology and the issue of the 
£-regime for staggered quarks is important because im- 
proved staggered fermions are currently the only lattice 
formulation that offers the prospect of a comprehensive 
study of the QCD axial anomaly free from the systematic 
effects of lattice spacing and unphysically large sea quark 
mass. 



Acknowledgments 

This research is part of the EU integrated infrastruc- 
ture initiative hadron physics project under contract 
number RII3-CT-2004-506078. 

We thank: Ph. de Forcrand for his topological charge 
measurement code; A. Hasenfratz for help in implement- 
ing the Hyp operator; CP. Lepage, P. Damgaard and 
G. Akemann for useful discussions. E.F., C.T.H.D. and 
Q.M. are supported by PPARC and A.H. by the U.K. 
Royal Society. The eigenvalue calculations were car- 
ried out on computer clusters at Scotgrid and the Dallas 
Southern Methodist University. We thank David Martin 
and Kent Hornbostel for assistance. 



[1] E. Follana, A. Hart and C. T. H. Davies, Phys. Rev. 

Lett. 93, 241601 (2004), [hep-lat/0406010]. 
[2] E. Follana, Nucl. Phys. Proc. Suppl. 140, 141 (2005), 

[hep-lat/0409062]. 
[3] K. Y. Wong and R. M. Woloshyn, hep-lat/0407003. 
[4] K. Y. Wong and R. M. Woloshyn, Phys. Rev. D71 

(2005), [hep-lat/0412001]. 
[5] M. Atiyah and I. Singer, Bull. Amer. Math. Soc. 69, 422 

(1963). 

[6] M. Atiyah and I. Singer, Ann. Math. 87, 596 (1968). 
[7] P. H. Damgaard, Nucl. Phys. Proc. Suppl. 106, 29 

(2002), [hep-lat/0110192]. 
[8] G. Akemann and P. H. Damgaard, Phys. Lett. B583, 

199 (2004), [hep-th/0311171]. 
[9] H. Leutwyler and A. Smilga, Phys. Rev. D46, 5607 

(1992). 



[10] E. V. Shuryak and J. J. M. Verbaarschot, Nucl. Phys. 

A560, 306 (1993), [hep-th/9212088]. 
[11] S. M. Nishigaki, P. H. Damgaard and T. Wettig, Phys. 

Rev. D58, 087704 (1998), [hep-th/9803007]. 
[12] P. H. Damgaard and S. M. Nishigaki, Phys. Rev. D63, 

045012 (2001), [hep-th/0006111], see corrections at hep- 

th/0006111v2. 

[13] J. J. M. Verbaarschot and I. Zahed, Phys. Rev. Lett. 70, 

3852 (1993), [hep-th/9303012]. 
[14] P. Forrester, Nucl. Phys. B 402, 709 (1993). 
[15] M. Goherman and J. Smit, Nucl. Phys. B245, 61 (1984). 
[16] J. Smit and J. C. Vink, Nucl. Phys. B286, 485 (1987). 
[17] S. J. Hands and M. Teper, Nucl. Phys. B347, 819 (1990). 
[18] L. Venkataraman and G. Kilcup, Nucl. Phys. Proc. 

Suppl. 63, 826 (1998), [hep-lat/9710086]. 
[19] A. Hasenfratz, "How good is the HYP staggered ac- 



17 



tion?" , Talk presented at Lattice 2003 (Tsukuba, Japan) . 

E. Follana et al, Nucl. Phys. Proc. Suppl. 129, 384 
(2004). 

M. E. Berbenni-Bitsch, S. Meyer, A. Schafer, J. J. M. 
Verbaarschot and T. Wettig, Phys. Rev. Lett. 80, 1146 
(1998), [hep-lat/9704018]. 

P. H. Damgaard, U. M. Heller and A. Krasnitz, Phys. 

Lett. B445, 366 (1999), [hep-lat/9810060]. 

M. Gockcler, H. Hehl, P. E. L. Rakow, A. Schafer 

and T. Wettig, Phys. Rev. D59, 094503 (1999), [hep- 

lat/9811018]. 

P. H. Damgaard, U. M. Heller, R. Niclasen and K. Rum- 
mukainen, Phys. Rev. D61, 014501 (2000), [hep- 
lat/9907019]. 

P. H. Damgaard, U. M. Heller, R. Niclasen and 
K. Rummukainen, Phys. Lett. B495, 263 (2000), [hep- 
Iat/0007041]. 

A. Hasenfratz, Phys. Rev. D64, 074503 (2001), [hep- 
lat/0104015]. 

MILC, C. Bernard et al, Nucl. Phys. Proc. Suppl. 119, 

991 (2003), [hep-lat/0209050]. 

MILC, C. Bernard et al., Phys. Rev. D68, 114501 (2003), 
[hep-lat/0308019]. 

MILC, C. Aubin et al., hep-lat/0409051. 

J. B. Kogut, D. K. Sinclair and M. Teper, Nucl. Phys. 

B348, 178 (1991). 

HEMCGC, K. M. Bitar et al, Phys. Rev. D44, 2090 
(1991). 

Y. Kuramashi, M. Fukugita, H. Mino, M. Okawa and 

A. Ukawa, Phys. Lett. B313, 425 (1993). 

B. AUes, M. D'Eha and A. Di Giacomo, Phys. Lett. 
B483, 139 (2000), [hop-lat/0004020]. 

S. Durr and C. Hoclbling, Phys. Rev. D69, 034503 

(2004) , [hep- lat/031 1002], 

S. Durr, C. Hoelbling and U. Wenger, Phys. Rev. D70, 

094502 (2004), [hep-lat/0406027]. 

S. Durr and C. Hoelbling, hep-lat/0408039. 

S. Durr, C. Hoelbling and U. Wenger, hep-lat/0409108. 

S. Durr and C. Hoelbling, Phys. Rev. D71, 054501 

(2005) , [hep-lat/0411022]. 

G. P. Lepage, Phys. Rev. D59, 074502 (1999), [hep- 
lat/9809157]. 

MILC, C. W. Bernard et al, Phys. Rev. D61, 111502 

(2000) , [hcp-lat/9912018]. 

MILC, K. Orginos, D. Toussaint and R. L. Sugar, Phys. 
Rev. D60, 054503 (1999), [hep-lat/9903032]. 
S. Naik, Nucl. Phys. B316, 238 (1989). 

F. Knechtli and A. Hasenfratz, Phys. Rev. D63, 114502 

(2001) , [hep-lat/0012022]. 

E. Follana et al., Further improvements to staggered 
quarks on the lattice, in preparation. 

G. Curci, P. Menotti and G. Paffuti, Phys.Lett. B130, 
205 (1983). 

Erratum:, G. Curci, P. Menotti and G. Paifuti, 
Phys.Lett. B135, 516 (1984). 

M. Luscher and P. Weisz, Commun. Math. Phys. 97, 59 

(1985). 

Erratum:, M. Luscher and P. Weisz, Commun. Math. 
Phys. 98, 433 (1985). 

M. Alford, W. Dirnrn, G. Lepage, G. Hockney and 
P. Mackenzie, Phys. Lett. B361, 87 (1995), [hep- 



lat/9507010]. 

[50] F. D. R. Bormet, D. B. Leinweber, A. G. Wilhams and 
J. M. Zanotti, Phys. Rev. D65, 114510 (2002), [hep- 
lat/0106023]. 

[51] J. B. Zhang et al., Phys. Rev. D65, 074510 (2002), [hep- 
lat/0111060]. 

[52] HPQCD, C. T. H. Davies et al, Phys. Rev. Lett. 92, 

022001 (2004), [hep-lat/0304004]. 
[53] P. de Forcrand, M. Garcia Perez and I.-O. Sta- 

matescu, Nucl. Phys. Proc. Suppl. 47, 777 (1996), [hep- 

lat/9509064]. 

[54] P. de Forcrand, M. Garcia Perez and I.-O. Stamatescu, 

Nucl. Phys. B499, 409 (1997), [hep-lat/9701012]. 
[55] J. CuUum and R. A. Willoughby, Lanczos Algorithms for 

Large Symmetric Eigenvalue Computations. Vol 1, theory 

(SIAM Society for Industrial & Applied Mathematics, 

Philadephia, 2002). 
[56] J. Hein, Q. Mason, G. Lepage and H. Trottier, Nucl. 

Phys. Proc. Suppl. 106 (2002), [hep- lat/01 10045]. 
[57] H. Trottier, G. Lepage, P. Mackenzie, Q. Mason and 

M. Nobes, Nucl. Phys. Proc. Suppl. 106 (2002), [hep- 

lat/0110147]. 

[58] W. Lee and S. Sharpe, Phys. Rev. D66, 114501 (2002), 

[hep-lat/0208018]. 
[59] C. Aubin et al, Phys. Rev. D70, 031504 (2004), [hep- 

lat/0405022]. 

[60] HPQCD and UKQCD, Q. Mason et al, Phys. Rev. Lett., 

in press , [hep-lat/0503005]. 
[61] L. Giusti, F. Rapuano, M. Talevi and A. Vladikas, Nucl. 

Phys. B538, 249 (1999), [hep-lat/9807014]. 
[62] UKQCD, D. A. Smith and M. J. Teper, Phys. Rev. D58, 

014505 (1998), [hcp-lat/9801008]. 
[63] R. G. Edwards, U. M. Heller, J. E. Kiskis and 

R. Narayanan, Phys. Rev. Lett. 82, 4188 (1999), [hep- 

th/9902117]. 

[64] R. G. Edwards, U. M. Heller, J. E. Kiskis and 
R. Narayanan, Phys. Rev. D61, 074504 (2000), [hep- 
lat/9910041]. 

[65] P. H. Damgaard, R. G. Edwards, U. M. Heller and 
R. Narayanan, Phys. Rev. D61, 094503 (2000), [hep- 
lat/9907016]. 

[66] P. Hasenfratz, S. Hauswirth, T. Jorg, F. Niedermayer 
and K. Holland, Nucl. Phys. B643, 280 (2002), [hep- 

lat/0205010]. 

[67] W. Bietenholz, K. Jansen and S. Shcheredin, ,]HEP 07, 

033 (2003), [hep-lat/0306022]. 
[68] L. Giusti, M. Luscher, P. Weisz and H. Wittig, JHEP 

11, 023 (2003), [hep-lat/0309189]. 
[69] P. H. Damgaard, private communication, 2004. 
[70] C. Aubin et al, hcp-lat/0402030. 

[71] B. Bunk, M. Delia Morte, K. Jansen and F. KnechtU, 
Nucl. Phys. B697, 343 (2004), [hep-lat/0403022]. 

[72] B. Bunk, M. Delia Morte, K. Jansen and F. Knechtli, 
hep-lat/0408048. 

[73] A. Hart and E. Miiller, Phys. Rev. D70, 057502 (2004), 
[hep-lat/0406030]. 

[74] D. H. Adams, hep-lat/0411030. 

[75] D. H. Adams, Nucl. Phys. Proc. Suppl. 140, 148 (2005), 

[hep-lat/0409013]. 
[76] F. Maresca and M. Peardon, hep-lat/0411029. 



